LAGUERRE FUNCTIONS AND REPRESENTATIONS OF 5u(l, 1) 



WOLTER GROENEVELT 

Dedicated to Tom Koornwinder on the occasion of his 60th birthday. 

Abstract. Spectral analysis of a certain doubly infinite Jacobi operator leads to orthogonality 
relations for confluent hypergeometric functions, which are called Laguerre functions. This dou- 
bly infinite Jacobi operator corresponds to the action of a parabolic element of the Lie algebra 
su(l, 1). The Clebsch-Gordan coefficients for the tensor product representation of a positive and 
a negative discrete series representation of su(f , f ) are determined for the parabolic bases. They 
turn out to be multiples of Jacobi functions. From the interpretation of Laguerre polynomials and 
functions as overlap coefficients, we obtain a product formula for the Laguerre polynomials, given 
by a discontinuous integral over Laguerre functions, Jacobi functions and continuous dual Hahn 
polynomials. 



1. Introduction 

Many special functions of hypergeometric type have an interpretation in representation theory of 
Lie groups and Lie algebras, see for example the book by Vilenkin and Klimyk [TTj- In this paper we 
consider the three dimensional Lie algebra su(l, 1), generated by H, B and C. Elements of su(l, 1) 
are either elliptic, parabolic or hyperbolic elements, which correspond to the three conjugacy classes 
of the Lie group SU {1,1). The self- adjoint element Xa = —aH + B — C, a G M, is an elliptic element 
for \a\ > 1, a parabolic element for \a\ = 1, and a hyperbolic element for \a\ < 1. In 11^ Koelink and 
Van der Jeugt consider the action ol Xa in tensor products of positive discrete series representations. 
This leads to convolution identities for several hypergeometric orthogonal polynomials. The idea 
to look at elements of the form of Xa is due to Granovskii and Zhedanov The action of Xa 
in the tensor product of a positive and a negative discrete series representation is considered in |7] 
for the elliptic case, and in for the hyperbolic case. In this paper we investigate the remaining, 
parabolic, case. In the quantum version of Xa is studied. The Lie algebra 5u(l, 1) is replaced by 
the quantized universal enveloping algebra Z^g(su(l, 1)) and Xa is replaced by a twisted primitive 
element. It turns out that in Z//g(su(l, 1)) the three cases are all the same. 

There are four classes of irreducible unitary representations of 5u(l, 1), the positive and negative 
discrete series, the principal unitary series and the complementary series. The tensor product of 
a positive and a negative discrete series representation decomposes into a direct integral over the 
principal unitary series. Discrete terms can occur, and these terms correspond to one complemen- 
tary series, or a finite number of discrete series. The Clebsch-Gordan coefficients for the standard 
bases are multiples of continuous dual Hahn polynomials. 

We consider the element X = —H + B — C, which is a parabolic element. In the discrete series 
X acts on the standard (elliptic) basis as a Jacobi operator, which corresponds to the three term 
recurrence relation for Laguerre polynomials. In the principal unitary series and the complementary 
series X acts on the standard basis as a doubly infinity Jacobi operator, which corresponds to the 
recurrence relation for Laguerre functions. So the Laguerre polynomials and functions appear as 
overlap coefficients between the (generalized) eigenvectors of X and the standard basis vectors. 
Using the differential equation for the Laguerre polynomials, we realize the generators H, B and C 
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in the discrete series as differential operators acting on polynomials. In these realizations the action 
of the Casimir operator can be identified with the hyper geometric differential equation, which leads 
to Jacobi functions as Clebsch-Gordan coefficients for parabolic basis vectors. This gives a product 
formula for Laguerre polynomials, which has a similar structure as the discontinuous integral for 
Bessel functions of Weber and Schafheitlin. 

This paper is organized as follows. In ^we consider a certain doubly infinity Jacobi operator, 
which corresponds to the action of X in the principle unitary series. Spectral analysis leads to 
orthogonality relations for Laguerre functions. This section is based on Jl] by Masson and Repka. 

In ^we turn to representations of the Lie algebra su(l, 1). We introduce the orthogonal poly- 
nomials and functions that we need in ^3.11 and we give some of their properties. In H3.2I we 
introduce the Lie algebra su(l, 1) and give the irreducible unitary representations. In ^3.31 we di- 
agonalize the element X in the various representations, and we give generalized eigenvectors. In 
^3.41 the generators B and C are realized as differential operators. Then the Casimir operator 
in the tensor product can be identified with the hypergeometric differential operator, and this gives 
Jacobi functions as Clebsch-Gordan coefficients. As a result we obtain a product formula for La- 
guerre polynomials, which involves Jacobi functions, Laguerre functions and continuous dual Hahn 
polynomials. 

Notations. If d^{x) is a positive measure, we use the notation dfj,2(^x) for the positive measure 
with the property 

(i(/U2 X fi^){x,x) = dfi{x). 
The hypergeometric series is defined by 

/ ai , . . . , Op \ _ (ai)n. . . . {ap)n 

where (a)„ denotes the Pochhammer symbol, defined by 

Y(a -I- n] 

(a)„ = ^— — ^ = a(a + l)(a + 2)...(a + n-l), nG 
r(a) 

For the confluent hypergeometric function we use the notation 

'a 



lFi a;6;z =1^1,, 
V 

and the second solution of the confluent hypergeometric differential equation is defined by 

^' ^) = w^^^ ill, 1-^1 («; ^) + ^-^i^z^-\F^{a _ 6 + l; 2 - 6; z), (LI) 
r(a — 0+1) r(a) 

see jl6l (1.3.1)]. This is a many valued functions of z, and we take as its principal branch that 
which lies in the complex plane cut along the negative real axis (— cxd, 0]. 

2. Laguerre functions 

In this section we determine the spectral measure of a certain doubly infinite Jacobi operator. 
This operator is obtained from the action of the self-adjoint element X of the Lie algebra 5u(l, 1) 
in the principal unitary series representation, see ^3.31 The eigenfunctions which are needed to 
describe the spectral measure, are called Laguerre functions. See |14j or ^U] for doubly infinite 
Jacobi operators. The calculation of the eigenfunctions and the Wronskian is obtained from (14j, 
but we repeat the calculations here briefly. 

The doubly infinite Jacobi operator L : ^^(Z) ^iJL) is defined by 

Lcfe = Ofeefc+i + ftfeefc + aA;_iefc_i, (2.1) 
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where {ek}kez is the standard orthonormal basis of £^(Z) and 

ak = ak{p,e) = \ {k + e + ^ + ip)\, 
bk = bk{p,e) = 2{k + e), 
where p>0,e£[0, 1) and (p,e) ^ (0, i). 

Remark 2.1. There exists a symmetry for the parameters of L. Let us denote L = L{p,£). The 
unitary operator U : Ck ^ (— l)'^e_fc intertwines L{p,e) with — L(— p, — e). So if f{z;e,p,k) is a 
solution to the eigenvalue equation Lf = zf, then {—\)^f{—z;—e,—p,—k) is another solution to 
the same eigenvalue equation. 

Observe that L is an unbounded, symmetric operator. The domain P of L is the dense subspace of 
finite linear combinations of the basis vectors e^. We define for a function / = X^^L-oo fk^k £ ^^(Z) 

oo 

L* f = ^ {akfk+i + bkfk + ak-ifk-i)ek- 

k=—oo 

on its domain 

V* = {f £ I L*f G e'^{z)} , 

then {L*,V*) is the adjoint of Note that = L. 

Solutions to Lv = —zv can be given in terms of confluent hypergeometric functions (see |lfij). 

Proposition 2.2. The following functions are solutions to Lv = —zv: 

|P('/j _|_ £ _|_ 1 _|_ 1 

Sk{z;p,e) = (-1)'=L-^— ^ m^Fi{k + £ + -+ip;l + 2ip;z), 

T{k + e + ^-tp) 2 

^ \T{^-k-e-ip)\ A , 
tk{z; p, e) = .^^ iFi{--k-e-ip;l- 2ip; -z), 

Uk{z;p,e) = {-l)''\r{k + e + ^ + ip)\U{k + e + ^ + ip;l + 2ip;z), z^ (-oo,0] 
Vk{z; p, e) = |r(^ -k-e-ip)\ - k - e - ip;l - 2ip; -z), z [0, oo). 

Proof. The first solution Sk follows from [TBI, (2.2.1)] 

(6 - a) iFi(a - l;b;z) + {2a - b + z) iFi{a;b;z) - aiFi{a + l;b; z) = 0. 

The second solution tk follows from the first using the symmetry relation for the parameters, 
cf. Remark 12.11 In the same way we find from 16, (2.2.8)] 

U{a - 1; b; z) - {2a - b + z)U{a; b; z) + a{a - b + l)U{a + 1; b; z) = 0, 

that Uk and Vk are solutions to Lv = —zv. □ 

The solution space to Lv = —zv is two-dimensional, since for a fixed n € v is completely 
determined by the initial values Vn-i and Vn- So the eigenfunctions given in Proposition 12.21 can 
be expanded in terms of each other. 

Proposition 2.3. We have the connection formulas 

Uk{z) =A{z)sk{z) + B{z)tk{z), z (-00, 0] 
Vk{z) =C{z)sk{z) + D{z)tk{z), z [0,oo), 
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where 



A{z) =T{-2ip), 

Biz) =.--W(2^p)-"""(^ + ^ + ^^) 



sm7r(e + — ip)\' 



C{z) ={-zf'Pe-'r{-2ip) 
D{z 

Or equivalently, for z £ C\M, 



sin7r(e + ^ + ip) 
sin 7r{e + ^ — ' 

D{z) =T{2ip). 



where 



Skiz) = E{z)uk{z) + F{z)vk{z), 
tk{z) = G{z)ukiz) + H{z)vk{z), 

E{z) = - sin^(e + \- ip)T{l + 2ip)e'^^^'+'^+'P\ 

TT 2 

F{z) = -- I sin7r(e + \ + ip)\ T{1 + 2ip)e^z-2ipgm«(e+i+ip)^ 
vr 2 

G{z) = - I sinvr(e + ^ + ip)\ T{1 - 2ip)e-^{-zf'Pe'''^'^'+^+^P\ 
vr 2 

H{z) = --sinvr(e + \- ip)T{l - 2ip)e'^^^'^^+'P\ 
vr 2 

where ^ = sgn9(z). 

Proof. The first connection formula follows from the reflection formula for the F-function, 

and Kummer's transformation: iF\{a;h]z) = e^iFi{b — a;b; —z). The second connection formula 
can be derived from the first using the symmetry for the parameters, see Remark 12.11 

The other two connection formulas follow from F16' (1.9.1), (1.4. 10)] or they can be derived from 
the first two. □ 

Definition 2.4. For two functions f{z) = Yl'kL-oo fk{z)ek and g{z) = Yl'kL-oo 9ki^)^k, the Wron- 
skian is defined by 

[f{z),g{z)]k = ak{fk{z)gk+i{z) - fk+i{z)gk{z)). 

If f[z) and g{z) are solutions to the eigenvalue equation Lv = zv, the Wronskian [f{z),g{zy\k is 
is independent of k, so the Wronskian can be found by taking the limit k — > ±00. Moreover, f{z) 
and g{z) are linearly independent solutions if and only if [f[z),g{z)] 7^ 0. 

Lemma 2.5. For < | arg(z)| < vr and k ^ 00 

Uk{z) = {-if ^ e^^-"/^^^' z--i-'Pk--i {l + 0(A:-i) 
v_k{z) = 0Fe-5^-V-4{fc-e>(_2)-i+ip^-i (^i + 

Sk{z) = -^^^ g|^+V4(fc+£)^^-3-vr(l + 2ip)k~\ (l + o{k-^2 
2vvr V 

t_k{z) = ^ g-^+V-4{fc-e)^(_^)-i+vr(i _ 2ip)k-'^ (l + 0{k-'2)^ . 
Proof. This follows from the asymptotic behaviour for |y| — > 00 of the modified Bessel functions 
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and the asymptotic expansions for the confluent hyper geometric functions in terms of modified 
Bessel functions [Ml (4.6.42), (4.6.43)]. □ 

For z £ C\R we introduce the spaces 

{oo oo 
fc=-oo k=0 ) 

S- = [f{z)= E fk{z)eu\L*f{z) = -zf{z)^nd E IA(^)|' < ooj . 

k fc=— oo k=—oo ) 

Then the deficiency space for L is H S~ . Note that dim < 2, and in case dim = 2, we 
have = S~ , since the sohition space to Lf = —zf is two dimensional. 

Next we put (pz = z^^u{z), z (oo,0], and = i—z)~^^v{z), z [0, cxd). From the transforma- 
tion U{a; b; z) = z^-^U{a - 6 + 1; 2 - 6; z) it follows that 



Ul{z) = Z^'Pui{z), Z0(-CX),O], 

Vk{z) = {-zf'p'^, z0[O,oo). ^'^ 



So we have {(l)z)k = {(t>z)k and {^z)k = {^z)k, and in particular {4>x)k G M for x > and {^x)k G 1^ 
for X < 0. Note that {(l)x)k and {^x)k are even in p. 

We calculate the Wronskian [cpz, ^z\- From Lemma 12.51 and = k + 0(1), for k oo, we find 
for z (— oo, 0] 

[Siz),u{z)]= hm -e'z-^-^'PT{l + 2ip)A;5 ( eV<k+e+l)z-^4(k+e)z _ ^^^ik+e)z-^^k+e+l)z\ _ 
k^oo 2 V / 

And since 



1 



gV4{fc+e+i)^-V4{'=+^)^ - eV4('=+^)^-V4{fc+^+i)-) = 2^1 (^1 + 0(A;-t) ) , A; ^ oo, 
we obtain 

[s{z),u{z)] = e^z-^'Pr{2ip + 1), z^ (-0O, 0]. 
Then we find from the connection formulas of Proposition 12.31 

[s{z),u{z)] = F{z)[v{z),u{z)], zeC\R, 

and this gives 

sin7r(e + ^ + ip)\ 
So we find that (pz and ^z are linearly independent. 

Proposition 2.6. For < | arg(z)| < vr, we have = span{<j)z} , and S~ = span{<l>2}, and L is 
essentially self-adjoint. 

Proof. From Lemma [2.51 we see that (pz = z^fu{z) G 5+ and $2 = {—z)~'^''v{z) G S~ , for < 
I arg(z)| < TT. Masson and Repka prove in jl4l Thm.2.1] that the deficiency indices of L are 
obtained by adding the deficiency indices of the two Jacobi operators obtained by restricting L 
to £^(Z>o) (setting a_i = 0) and to £^(— N) (setting ao = 0). Since J2T=o^/^k and J2k=-oo^/'^k 
are divergent, -S, Ch.VII, Thm.1.3] proves that have deficiency indices (0,0), and hence so has 
L. So dim = 1, and the proposition follows. □ 



z, ^z] = z''{-zrP[u{z),viz)] = < I arg(z)| < vr. (2.3) 
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We use the Stieltjes- Perron inversion formula, see |^ §XII.4], to calculate the spectral measure; 
EfJia,b)) =limlim-^ T \g{x + i5)f,g) - {G{x - i5)f,g)dx. 
In this case the resolvent G{z) can be calculated explicitly by 

{G{z)f, g) = — + - \^ik). (2.4) 



Proposition 2.7. The spectral measure for the operator —L defined by (|2.1() . is described by the 
following integral, for B '^M a Borel set, 

{E{B)f,g) =\ [ e-lsmnie + ^ + ip)\\f,v{x)){v{x),g)dx 



/sn(-oo,o) 

\ {' 1 

+— / e"''|sin7r(e + - + z/9)p(/,u(x))(u(x),ff)cix, 

/Bn(o,oo) ^ 



^2 



where f,g £ £^(Z) and with the notation of Proposition[ 

oo oo 

u{x) = ^ Uk{x)ek, v{x) = ^ Vkix)ek. 

k=—c>o k=—oo 

Moreover, is not contained in the point spectrum of L. 
Proof. We define 



A, X ,. i^x+i5)k{4>x+i5)l {^x-i5)k(,(t'x-i5)l 

A(a;j = hm — — , (2.5) 

510 [(I):r+i5,^x+i5\ [4>x-i5,^x-i5\ 

then we have, using (|2.3() and S, = sgn9(2;), 

A(x) = - lim - I sin^(e + l + ip)\ (e'''^'+^\^^+^s)k{4>x+t5)i - e-'''^'+^\'^^-i5)k{(t>x-i5)i) ■ 

510 TT Z \ / 

For X < we have limsio{^x+i5)k = ^^''^5io{^x-i5)k = {^x)k, and from the connection formulas we 
find 

hm e^-(^+5)(0,+,,), - e--(^+5)(0,_.,), 

o|0 



lim fe^'^(^+l)(x + i5)^M(x + i(5)s;(3; + i6) - e-'^'^'+^^x - iSf^Aix - i6)si{x - i6)) 
5iO V / 

+ lim fe*^(^+5)(x + i6yPB{x + i6)ti{x + i6) - e-*^(^+5)(x - i6yPB{x - i6)ti{x - i6)) 

510 V / 



<5iO 
li 

510 

1 1 

■ 2isimT{e + - + ip){-xypr{-2ip)si{x) + 2ie''{-x)-'Pr{2ip) | sin7r(e + - + ip)| t/(x) 

: 2ie''(-x)~*^| sin7r(e + ^ + ^p)| ■y«(^)- 



Here we used 



hm(-y ± i(5)" = e±'™y", y > 0. 



For X > we use the symmetry for the parameters, cf. Remark 12.11 So we find 

2i 1 

— e^(-x)"2*^| sin7r(e + - + ip)\^Vk{x)vi{x), x <0, 

TT 2 



A(x) = < 



(2.6) 



-— e '^x^*''! sin7r(e + - + ip)ptifc(x)u;(x), x > 0. 
vr 2 
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Both expressions are clearly symmetric in k and Z, so the sum in (|2.4)) can be antisymmetrized 
using (|2.2j) . Now, if is not contained in the point spectrum of L, the result follows from the 
Stieltjes- Perron inversion formula and (|2.2|) . 

To show that is not an element of the point spectrum of L, we show that kerL = {0}. First we 
calculate the Wronskian [s(0),t(0)], using Definition 12.41 with A; = 0. A straightforward calculation 
gives 



[s{Q),m] = 2ip 



sin7r(e + ^ — ip) 
I sin7r(e + i + ' 



hence s(0) and t(0) are linearly independent. So if / G kerL, f ^ 0, then / is a linear combination of 
s(0) andt(O). But s(0),t(0) since |sfc(0)| = 1 and \tk{0)\ = 1, and therefore kerL = {0}. □ 

Remark 2.8. The result of Proposition lTTl remains valid if is replaced by where A G (~|) 
and e G [0, or A G (— ^, e — 1) and e G 1). In this case operator L is obtained from the action 
of the self-adjoint element X in the complementary series representation of 6u(l, 1), see ^3.31 

Let us define the Laguerre functions ipnix;p,£), n G Z, by 

Tpn{x;p,e) -- 



Vn{x;p,£), X < 0, 
Un{x;p,e), X > 0, 



and we define the weight function w{x; p, e) by 

w{x;p,e) = ^1 sin7r(e + ^ + i/))pe~l^l. 

TT^ 2 
From Proposition 12 . 71 we find the following. 

Theorem 2.9. For p > 0, e G [0, 1) and (p, e) (0, the Laguerre functions il^n{x] p, e) form an 
orthonormal basis of Lp'{^,w{x; p,e)dx). 

Proof. The orthonormality of the Laguerre functions follows from Proposition 12.71 by replacing / 
and g by standard orthonormal basis vectors e„ and Cm, n,m £ Z, and using ;B = M. Completeness 
of the Laguerre functions follows from the uniqueness of the spectral measure. □ 

We also want to define the Laguerre functions for x = 0. In order to find the natural definition 
in X = we calculate A(0), where A(x) is defined by (|2.5|) . Using the connection coefficients from 
Proposition 12.31 we find 

A(0) = -hm i |sin7r(e + i + zp)|(e*-(^+5)($^,),(</>,,), _ e-'-^'+'^\<^^,s)k{4>~^s)l 
510 TT 2 V 

= - lim -\ sin^(e + - + ip)| f e^'^^^+s+^p) \C{i6)sk{iS) + D{i5)tk{iS)] \A{i5)si{i6) + B{i5)ti{i5)] 

SiO TT 2 \ 



e 



-Me^+ip)^Ci-i6)sk{-i6) + D{-i6)tk{-iS)] [A{-i6)si{-i6) + B{-i6)tii-i6)]y 
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To compute this limit, we use 

limskiiS) = limski-id) = Sfc(O), 
sio 510 

limtkiiS) = \imtki-i5) = tfc(O), 
limA{i5)D{i6) = limA{-i6)D{-i6) = \T{2ip)f , 

510 sin7r(e + ^ — ip) 



lim B{-i6)C{-i6) = e'^'^'P \T{2ip)\ 



2 sin7r(e + 2+ ip) 



510 sin7r(e + ^ — ip) ' 

then we find 

A(0) = sin^(e + ^ + ip) I sin7r(e + ]^ + ip)\ \T{2ip)\^ (tfe(O)sKO) + SkWM) ■ 
From Euler's reflection formula for the F-function we obtain 



, , sinvrfe + i — zp) — -— 
I sm7r(e + 2+ip)\ 



and this gives 



A(0) = -- I sin^(e + i^ + ip) T{2ip)\'' {tk{f))tm + tfc(O)tKO)) 



Comparing this result with (|2.6j) . we see that for x = the Laguerre function can be defined by 

^/;„(0;p,e) = |F(2ip)| 



tk{o)J • 



3. Clebsch-Gordan coefficients for parabolic basis vectors of su(1, 1) 

3.1. Orthogonal polynomials and functions. The Wilson polynomials, see Wilson jl9j or ^ 
§3.8], are polynomials on top of the Askey-scheme of hypergeometric polynomials, see Koekoek 
and Swarttouw The continuous dual Hahn polynomials are a three-parameter subclass of the 
Wilson polynomials, and are defined by 

Sn{y;a,b,c) = {a + b)n{a + c)n3F2\ . , ; 1 , x = y. (3.1) 

\ a + b,a + c J 

For real parameters a, 6, c, with a + 6, a + c, b + c positive, the continuous dual Hahn polynomials 
are orthogonal with respect to a positive measure, supported on a subset of M. The orthonormal 
continuous dual Hahn polynomials are defined by 

Q ( , N (-l)"sn(?/;a,^,c) 



yjn\{a + b)n{a + c)„(6 + c)„ 



By Kummer's transformation, see e.g. 01, Cor. 3.3.5], the polynomials Sn and Sn are symmetric in 
o, b and c. Without loss of generality we assume that a is the smallest of the real parameters a, b 
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and c. Let a, b, c) be the measure defined by 



/ 

JR 



r{b-a)r{c-a) 
^ r(-2a)r(6 + c) 



r(a + ix)T(b + ix)T{c + ia;) 



r(2ix) 



fix' 



T{a + b)T{a + c)T{b + c] 



dx 



j^Nfc (2Q)fc(Q+ l)fc(Q + b)fc(a + c)fc 
f-^ ^ (a)fc(a-6+l)fc(a-c + l)fcA;!-^^ 



-ia + kf), 



where K is the largest non- negative integer such that a + K < 0. In particular, the measure 
dfi{-; a, b, c) is absolutely continuous if a > 0. The measure is positive under the conditions a+b > 0, 
a + c > and 5 + c > 0. Then the polynomials Sn{y', a, b, c) are orthonormal with respect to the 
measure dfi{y; a, b, c). 



The Laguerre polynomials are defined by 

(a + l)r 



The orthonormal Laguerre polynomials 



n! 



iFi{-n; a + l;x) 



(3.2) 



n! 



V (« + l)n " 

are orthonormal on [0, oo) with respect to the weight function 

^ r(a + i)- 

They satisfy the three-term recurrence relation 

xli^\x) = -V(n + l)(a + n+l)4'5i(a;) + (2n + a + 1)4°^^) " V^in + a)i%ix), 
and the differential equation 

xy"{x) + (« + ! — x)y'{x) + ny{x) = 0, y(x) = /^°^(x). 
The Jacobi functions, see J2]) are defined by 



^\ '{x) = 2F1 



i(a + /? + 1 - iX), i(a + /3 + 1 + a) 
a + 1 



(3.3) 



Here we use the unique analytic continuation to C \ [l,oo) of the hypergeometric function. The 
Jacobi functions are eigenfunctions of the hypergeometric differential operator 

d^ d 
-x(l + x)— J - [a + 1 + (a + /? + 2)3;] — 
ax'^ dx 

for eigenvalue \[{a + (3 -\- Vf' + A^]. Spectral analysis of the hypergeometric differential operator 
leads to a unitary integral transform called the Jacobi-function transform. The Jacobi-function 
transform is given by 



{Tf){\) = / /(x)<^^'')(x)A„,^(x)da 



(3.4) 
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where a > —1, /3 G M, Aa,i3{x) = 2^"+^^+^x"(l + x)^, and di'{X) is the measure given by 
1 /" 1 

— / g{X)di'{X) = — / g{X)\ca^/3{X)\-'^dX-i^g{X)Res{ca,|3{^J■)ca,|3{-^J■))~\ 

J Jo A* 

2°+^+i-^^r(a + i)r(a) 



r(i(a + /3 + 1 + iA))r(i(a - /? + 1 + iX)) ' 
P = - a - 1 - 2j) I j G Z>o, - a - 1 - 2j > o}. 

Observe that the measure dv{X) is absolutely continuous if < a + 1. 



3.2. The Lie algebra su(l, 1). The Lie algebra su(l, 1) is a three dimensional Lie algebra, gen- 
erated by H, B and C satisfying the commutation relations 

[H, B] = 2B, [H, C] = -2C, [B, C] = H. (3.5) 

There is a *-structure by H* = H and B* = —C. The Casimir operator J7 is a central element of 
U {su{l, 1)) , and Q is given by 

n = -^{H^ + 2H + 4CB). (3.6) 

There are four classes of irreducible unitary representations of 5u(l, 1), see §6.4]: 
The positive discrete series representations vr^ are representations labelled by A; > 0. The 
representation space is ^^(Z>o) with orthonormal basis {e„}„gz>o- The action is given by 

TT^iH) e„ = 2(/c + n) e^, 
7r+ {B)en = V(n + l)(2/c + n) e„+i , 
7r^(C)ert= - yjn{2k + n - 1) e„_i, 
7r^(ri) e„ = /c(l - A;) e^. 

The negative discrete series representations vr^ are labelled by > 0. The representation space 
is £^(Z>o) with orthonormal basis {e„}„gz^,j. The action is given by 

'K~^{H)en= -2{k + n)en, 

TT7{B)en= - yjn{2k + n - 1) e„_i, 

_ (3.8) 

^fc (C) en = \/(n + l)(2/c + n) e^+i, 

7r^(ri) e„ = k{l - k) en- 

The principal series representations vr^'^ are labelled by e G [0, 1) and /O > 0, where (p, e) ^ (0, ^). 
The representation space is ^^(Z) with orthonormal basis {e^lnez- The action is given by 

■kP^'{H) e„ = 2(e + n)en, 

7r^'^(5) e„ = |(n + e + i + ip)| e„+i, 

1 (3.9) 

7r^'^(C)e„= -|(n + e-- + ip)|e„_i, 

7r'''^(0)e„= (p2 + i)en. 
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For e) = (0,^) the representation vr'^'a splits into a direct sum of a positive and a negative 
discrete series representation: vr'^'a = 7r|' © vr^. The representation space splits into two invariant 

2 2 

subspaces: {e^ | n < 0} © {e„ | n > 0}. 

The complementary series representations vr^'^ are labelled by e and A, where e E [0, ^) and 
A G (— 1> — e) or e G (i, 1) and A S (— — 1). The representation space is ^^(Z) with orthonormal 
basis {en}n&- The action is given by 

TT^'^{H) e„ = 2(e + n) e^, 

^^'"(B) e„ = V(n + e + 1 + A)(n + e - A) e^+i, 

(3.10) 

^^•^(C)e„= - V(n + e + A)(n + e - A - 1) e„_i, 

7r^'^(0)e„= -A(l + A)en. 

Note that formally for A = —i^ + ip the actions in the principal series and in the complementary 
series are the same. 

We remark that the operators ()3.7() - (|3.10|) are unbounded, with domain the set of finite lin- 
ear combinations of the basis vectors. The representations are *-representations in the sense of 
Schmiidgen [El Ch.8]. 

The decomposition of the tensor product of a positive and a negative discrete series representation 
of su(l, 1) is determined in [3 Thm.2.2], see also XI , §8.7.7] for the group SU{1, 1). 

Theorem 3.1. For ki < k2 the decomposition of the tensor product of positive and negative discrete 
series representations o/su(l,l) is 

oo 

f® 11 

< ®^k,- I ^^''<^P^ h-k2> --, ki+k2>-, 





oo 



fc^wMfc^- , 7r'''^d/3 © 7r^'% ki + k2<^, 





oo 



1 1 



iez>o 

where e = ki — k2 + L, L is the unique integer such that e € [0, 1), and A = —ki — ^2- Further, 
under the identification above, 



Cm e„2 = (-1)"^ / Sn{y;ni - n2)em^^n2-Ldp'2{y;ni - n2), (3.11) 
JR 

where n = min{ni,n2}, Sn{y',p) is an orthonormal continuous dual Hahn polynomial, 

\ Sn{y; ki - k2 + ^,ki + k2 - ^,k2 - ki - p + i), p<0, 
\Sniy;k2 - ki + ^,ki + k2- l,ki-k2+p+^), p> 0, 



Sn{y;p) 

and 

dp{y;p) 



dp{y; ki - k2 + ^, ki + k2 - I, k2 - ki - p + i), p <0, 
dfi{y; k2 - ki + ^, ki + k2 - ^, ki - k2 + p + ^), p>0. 
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The inversion of (|I-{.1H1 can be given explicitly, e.g. for an element 



/ (g) Er-L 



f{x)er-Ldx e L'^iO,oo)<E)f{Z) 



f{Z)dx 



in the representation space of the direct integral representation, we have 

/ oo r „ 

oo p „ 

Y^i-ir / Sp(y;r)f{y)df,-2{y;r) 



/ (g) Cr-L = ' 



r < 0, 
r > 0. 



(3.12) 



For the discrete components in Theorem 13.11 we can replace / by a Dirac delta function at the 
appropriate points of the discrete mass of dfi{-;r). In the following subsections we assume that 
discrete terms do not occur in the tensor product decomposition. From the calculations it is clear 
how to extend the results to the general case. 



3.3. Parabolic basisvectors. We consider the self-adjoint element 

X = -H + B-C e5u{l,l), 

which is a parabolic element. We determine the spectral decomposition of X in the various repre- 
sentations. We also give (generalized) eigenvectors of X. This is done in the same way as in jllj . 
using (doubly infinite) Jacobi operators. First we consider X in the discrete series. The action of 
X can be identified with the three-term recurrence relation for the Laguerre polynomials. 

Proposition 3.2. The operators 0^ defined by 

: e\Z>o) ^ LmO,oo),w'^^''-'^\x)dx) 

are unitary and intertwine tt^{X) with Mzpx- 

Here M denotes the multiplication operator: Mfg{x) = f{x)g{x). 
In terms of generalized eigenvectors, Proposition 13.21 states that 

oo 
n=0 

is a generalized eigenvector of 7r^(X) for eigenvalue ^x. These eigenvectors can be considered as 
parabolic basis vectors for su(l, 1). 

Next we consider X in the principal unitary series. We find that ■kP'^[X) extends to a doubly 
infinite Jacobi operator which corresponds to the recurrence relation for the Laguerre functions. 
The spectral analysis of 7r'''^(X) is carried out in ^ 

Proposition 3.3. The operator QP''^ defined by 

eP'^ : f{Z) L^{R,w{x;p,e)dx) 

e„ ^ (-l)"V'n(-;P,e) 

is unitary and intertwines ttP'^{X) with M^- 
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So, for X E M, 

oo 

vP^^{x)= J2 {-irMx;p,e)en, 

n=~OD 

is a generalized eigenvector of itP''^{X) for eigenvalue x. Note that v'''^{0) is two dimensional, since 
tpn{0) is given by a two dimensional vector. 

Next we consider the action of X in the tensor product. Recall that A(Y) = 1 <^Y + Y fSil for 
Y £ su(l, 1). Then we find from Proposition 13.21 the following. 

Proposition 3.4. The operator T defined by 

T : f{Z>o)(E)f{Z>o) ^^([0,oo) X [0,oo),w^^'''-^\xi)w^^''^-^\x2)dxidx2) 

is unitary and intertwines vr^^ ® 7r^^(A(X)) with M^j-xi- 
So 

oo 

v+{xi) v~{x2) = Yl '^'^(2;2) en, e,, 

ni, 712=0 

is a generalized eigenvector of vr^^ 7r^^(A(X)) for eigenvalue X2 — xi. 

3.4. Clebsch-Gordan coefficients. We want to determine the Clebsch-Gordan coefficients be- 
tween the uncoupled eigenvectors v~^{xi)(^v~{x2) and the coupled eigenvectors /® v^^'^^xi — X2)dp. 
This comes down to finding the function for which the operator Tg defined by 

oo oo 

: l2(0,oo) ®£2(Z) ^ / f{Z)dx^ / {R,w{t; p,e)dt) dp, 



/■oo 

/ e„ (-1)" / f{p)g{p)'ipn{t] p, e)dp, 
Jo 

is the same as T. From Theorem 13.11 and (|3.9|) we know that the Clebsch-Gordan coefficient g must 
be an eigenfunction of the Casimir operator 17 in the tensor product for eigenvalue p"^ + j. So first 
we determine the actions of the generators H, B and C on parabolic basis vectors. 

We start with a very simple lemma, which is based on the fact that sl(2, C) is semi-simple, so 
[s[(2,C),sl(2,C)] =s[(2,C). 

Lemma 3.5. 

B = -\H,X] + -X + -H, C = -\H,X]--X --H. 
4^ ' ^ 2 2 ' 4^ ' ^ 2 2 

Proof. From the definition of X and the commutation relations (|3.5|) we find [H, X] = 2B + 2C. 

This proves the lemma. □ 

This lemma shows that to find the action of the generators H, B and C, it is enough to find the 
action of H, since the action of X is known. 

Proposition 3.6. In the positive discrete series, the generators H , B, C have a realization as 
differential operators acting on polynomials: 

d^ d 
7r+{H) = -2x—. - 2{2k - x)— + 2k, 
ax^ dx 



d^ 2(fc ^ 
dx"^ dx 
d^ . d 



TT+iB) = -X— - 2{k - x)— + {2k - x), 
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In the negative discrete series H , B, C have a realization as differential operators on polynomials: 



vr-(C) = -x^-2(fc-x)A + (2fe-^). 

Proof. We show that 0^ intertwines the actions of H, B and C given by ()3.7() . with the differential 
operators given in the proposition. 

From the differential equation for the Laguerre polynomials, we find for the action of H 

Q+7T+{H)en = {2n + 2k)l^^-'-\x) 

= -2x^e-'\x) - 2{2k - x)^e-'\x) + 2ke-'\x) 

- 2x^ - 2{2k -x)-^ + 2k) e+e„. 
ax^ dx J 

So we have realized it^(H) as a differential operator. By Proposition 13.21 Trt (X) is realized as the 
multiplication operator M-x- A direct calculation shows that 

e+TT+{[H,X])en = (4x^ + 2(2A; - x))e+e„. 

Then Lemma 13.51 proves the proposition for the positive discrete series. We find the action in the 
negative discrete series in the same way, or we use the Lie-algebra isomorphism i9, given by 

'd{H) = -H, ^{B) = C, ^{C) = B. 

Then 7r+(t?(y)) = tt^{Y) for Y G su(l, 1). 

A straightforward calculation shows that these operators indeed satisfy the 5u(l, 1) commutation 
relations. □ 

It is also possible to find the actions of H, B and C on the eigenvectors v''''^{x). This is done 
using the differential equation for the Laguerre functions, which follows from the confluent hyper- 
geometric differential equation. We do not need these actions here. 

Next we want to calculate tt^^ (8) vr^^ (A(il)) for these realizations. From (|3.(i|) we obtain 

A(J]) = 1®9. + 9.®1-]^H®H -{C®B + B®C). (3.13) 

Proposition 13. HI shows that vr^^ vr^^ (A(0)) is a differential operator acting on polynomials in two 
variables. Let p{xi) and q{x2) be polynomials in xi, respectively X2, then we find after a long 
calculation in which many terms cancel, 

-xiX2(p"{xi)q{x2) + 2p'{xi)q'{x2) + p{xi)q" {X2)) 

^ \ (3.14) 

+(2xiX2 - 2^1X2 - 2A;2Xi)(p'(xi)q'(x2) + p{xi)q' {X2) 



+ [2kiX2 + 2k2Xi - 2kik2 - X1X2 + ki{l - ki) + A;2(l - k2))p{xi)q{x2), 

where p'{xi) = -;^p{xi) and q'{x2) = ■^q{x2)- We show that vr^^ vr^^ (A(ri)) can be identified 
with the hyper geometric differential operator, and therefore the eigenfunctions are Jacobi functions. 
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Let i G R and define 



P{x,t) 

Proposition 3.7. The operator H, defined by 

E : P{x,t) H 



p{x — t)q{x), t < 0, 
p{x)q{x), t = 0, 
^p{x)q{x + t), t>0. 



(2fc2-l,2fci-l) 



i-x/t), t < 0, 
t = 0, 
t > 0, 



intertwines vr^^ VTj^.^ (A(ri)) with Mpi^ij^. 

Proof. First we assume t > 0. Put x = t = X2 — xi and e^ip{x] t) = p{x)q{x + 1) in H3.14() . tlien 

f - + t)^ - {2ki{x + t) + 2k2x)^ + - h) + A:2(l - fca) - 2A:iA:2)v? . 

L ax^ ax J 

Comparing this to the differential equation for the Jacobi functions, we see that the Jacobi function 
(f'^p^ ^'^'^^ ^^(^/^) is an eigenfunction of this differential operator for eigenvalue + |. And then 
the intertwining property of H follows. The case t < is proved similarly. 

Next we assume t = 0. We put x = xi = X2 and e^ip{x) = p{x)q{x) in 1)3. 14() . and then from 
Theorem 13. II it follows that we must solve the following Euler differential equation: 



{2ki + 2A;2)xy' + {ki + /c2)(l - h - k2)y = {p^ + ^)y. 



The general solutions to this equation are given by y = cix 



□ 



For t / we can identify the spectrum of vr^^ (8> 'k^_^[/S.{^)) with the support of the measure 
di'{2p). Note that, naturally, the support is exactly the same as the support of the orthonormality 
measure for the continuous dual Hahn polynomials given in Theorem 13.11 

Since the Jacobi function transform is unitary, there exists a constant c such that c H is unitary for 
t 7^ 0. So Proposition 13.71 determines for t ^ Q the Clebsch-Gordan coefficients for the eigenvectors 
up to a factor independent of x. To determine the exact Clebsch-Gordan decomposition for the 
eigenvectors, we need to find this factor. For t = the Laguerre functions are C^-valued, so the 
Clebsch-Gordan coefficients for t = are also C^-valued. We show that in this case the unitarity 
of c H corresponds to the unitarity of the Mellin transform. 

Theorem 3.8. The Clebsch-Gordan coefficients for the parabolic bases are given by 



?(xi - X2) 



(2A;2-l,2fci-l) 



X2 



Xi - X2 



-ki-k2+ip^ 



^ —ki-k2-ip 



Coip) e-^xl 
Mp) e^'xl 



n t \ ^1 (2fci-i,2fc2-i) 
C+{p)e V2p 



Xl 



X2 - Xl 



X2 - Xl < 0, 
X2 - Xl = 0, 
X2 — Xl > 0, 
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where 



C-ip) 
Coip) 

c+{p) 



27r 



lT{2ki 



V{ki + k2-\ + ip) 



T{ki -k2-\ + ip)T{2ip) 



_1 (-1)^ v/2r(2fci)r(2fc2) r(fci + fca - ^ + ip) 

|r(A:i + /C2 - ^ + ip)^{ki - /C2 + i + ip)\T{k2 -ki + ^-ip)' 

T{ki + k2-^ + ip) 



— (x2-2;i)'=^+'=2-5-^^ 



27r 



/r(2fc2 

T{2ki 



T{k2 -ki + ^+ ip)T{2ip) 



Proof. Recall that the Clebsch-Gordan coefficients are the functions g such that T = T^. Put 
X = xi, t = X2 — xi and assume t > 0. Since the Clebsch-Gordan coefficients do not depend on ni 
and n2 it is enough to show that for the function g defined in the theorem, we have T(eo (SD cq) = 
Tg ( Jq°° e-Ldp}/"^^ , where dp is the orthogonality measure for the continuous dual Hahn polynomials 
as in Theorem 13.11 Explicitly, we must prove the following identity 



1 



-11 



C+{p) 



^T{2ki)T{2k2 



T{ki + k2-^ + ip)T{k2 - ki + ^ + ip)T{ki -k2 + ^ + ip) 



T{2ip) 



dp. 



We use an integral representation for the second solution of the confluent hypergeometric differ- 
ential equation, see ^1 (3.2.55)], 



^a— c /"oo 



U{a] b; z) 



-zy c—l 



r(c) Jo 

with parameters given by 



a,l + a 
c 



a = ki - k2 + - + ip, 



b = l + 2ip, 



■,-y]dy, 3f?c> 0, JRz > 0, 



z = t. 



2ki, y 



By the definition of the Laguerre function V'n(i) for t > 0, see fjH the definition of a Jacobi function 
1)3. 3|) and Euler's transformation '1, (2.2.7)], we have 

ip-L{t;p,ki -k2 + L) = 
(-1)^2 



Lr)'i—Aki—ik2^\-k\-k2+ip 



jT{ ki-k2 + \+ip)\ 

T{2ki) 

Taking the inverse Jacobi transform of this, gives 



X' (2fci-l,2fe2-l) /a: 



- ) A2fci-1,2A;2-1 



x\ dx 



tJ t 



1 = (-1) 



27r 



^.^(2^,-1,2.2-1) t-^f>^_^{t- p, k,-k2+ L) 



\T{ki -k2 + ^ + ip)\ 
f(2h) 



r(A;i + k2-^ + ip) 



T{2ip) 



dp 



This is the desired identity. For t = xi — X2 < the theorem is proved similarly. 

For X — xi — X2 we determine the (C^-valued) function g for which we have Ti^em 
T3((-l)"2 Snem-n2-Ldp^^'^), see Theorem KT\ We start with HEl (3.7.4)]; 



-a;' s— 1 



-Fi(o; b; x)iFi{c; d; x)dx 



T{s)T{d)T{d-c-s) 
T{d- c)T{d- s) 



iF2 



a, .3,1 + s — d 
bA + s + c — d 
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where K(s) > and we assume a to be a non-positive integer. We consider this as a Melhn 
transform, and then the inverse transform gives, for A; > 0, 

1 /■^+'°° r(.)r(d)r(d-c-.) ( g^sA + s-d \_ _ 

— rprri ^ 3i^2 , ^ , , ]l]x ds = e iFi{a;b;x)iFi{c;d;x). 

2v^^ A-ioo T{d- c)T{d- s) \b,l + s + c- d J 



We use 



-ni, b = 2ki, c = -n2, d = 2k2, s = ki + k2 - ^ + ip, 



and then we transform the terminating 3i^2-series using the first formula on page 142 in Next 
we add to the integral the same integral with p replaced by —p, then we obtain an integral for which 
the integrand is even in p. Using Euler's reflection formula and the definitions for the Laguerre 
polynomials (|3.2|) . the continuous dual Hahn polynomials (|3.1I) . and the Laguerre functions of 
argument 0, we have 

1 f°° 11 1 

- / Sni{p'^;ki - k2 + -,ki + k2 - -,k2 - ki +n2 - ni + -) 
TT Jq 2 2 2 

^ /^n, ^ , l^(^2 - fci + :t + n2 - ni + 

\/ {2k2)n2-nA"^2 - nx)\ ^ 

X { ^^^^^^±-M±^t„,_„,_,(0; p, kx-k2^ L) x^^-^-- 



+ 



r(/c2 - + \ - ip) 

T{ki +k2-\-ip) 
T{k2 -ki + \ + ip) 



ini-n2-L(0; p, ki-k2+ L) x^^'^^-'^^+V 



Writing the last integral as 

(-1)"' 



gives the expression for the Clebsch-Gordan coefficients. 



AO) 



□ 



Remark 3.9. The explicit expressions for the Clebsch-Gordan coefficients as 2-^1-series can also be 
found in Basu and Wolf pj. The method used in |2j to compute the Clebsch-Gordan coefficients is 
different from the method used here. 



Theorem 13.81 gives the following Clebsch-Gordan decomposition for the parabolic basis vectors 



v'^{xi)^v {X2) 







X2 



Xi- X2 



V^'^{X2 — Xi)dp, X2 — Xi < 0, 



Coip) e^-xl 
\Co{p)e''x\ 



-ki-k2-ip 



vP'%0)dp, 



X2 - Xl = 0, 



Xl 



10 ■ \x2-x1 

From Theorem 13 . 81 we obtain a product formula for Laguerre polynomials. 



V^'^{X2 — Xl)dp, X2 — Xl > 0, 
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Theorem 3.10. The Laguerre polynomials satisfy the following product formula: 
for xi > X2 



1 r°° 1 1 



^ X2 (2fc2-l,2fci-l) / X2 



Xi - X2 



2^k2-h + -, 



U{n2 - ni + k2 - ki + - - ip;l - 2ip; xi - X2)dp, 



for xi < X2 



1 111 

4f l-l)(xi)42'=2-l)(x2) = — / d+ Sndp^; k2-k,+n2-ni + -, k^ + k2 - -, k^ - k2 + -, 



2tt 



X e^^ifi^^ ^''^^'^ { — — — ]U{ni-n2 + ki-k2 + ^+ip;l + 2ip; X2 - xi)dp, 
\X2 — xi ' ^ 



for xi = X2 = X 

1 111 

4f i-l)(x)Li2'=^-^)(x) = - / Sn.ip^; k2-k,+n2-m + -,k^+k2--,k,-k2 + - 

TT Jq III 



where 



d- 
do 



(xi - X2) 



ni\n2\T{2k2) 

{X2 - 

ni\n2\T{2ki) 



+ dox^-^^~'"'+'P 



T{ki +k2-\+ ip)T{k2 - ki + n2 - ni + ^ + ip) 



dp, 



T{2ip) 

T{ki + k2 - 1+ ip)T{ki - k2 + ni - n2 + \ + ip) 



T{2ip) 

11 1 

T{ki + k2 - - + ip){k2 - ki + - - ip)n2-ni- 



Proof. This follows from writing out explicitly T(e„^ CSCnj) = Tg( J® Sn{p'-,ni — n2)eni-n2-Ldp'2) , 
n = miii{ni, 71,2}, where g is given in Theorem 13.81 and using [3 (3.13)] 



1 



1 



Snip ]ki - k2 + -,ki + k2 - -,k2 - ki - p + 



1. 



(-lf\{ki - /c2 + ^ + ip)p\^Sn-p{p^]k2 - A;i + + A:2 - ^,^1 - ^2 + P+ ^) 



□ 



Remark 3.11. (i) The confluent hypergeometric ^/-function can be considered as a Whittaker func- 
tion of the second kind, see ^1 (1.9.6)]; 



11 1 

Wk,rn{x) = e"5^x'^+^C/(m - A: + -; 2m + 1; x) 



These Whittaker functions are the kernel in the Whittaker function transform, given by 



(W/)(A) 



27r 



f{x)Wk,ix{x)x 2dx, 

\wf){\)x-'2WkM^) 



r{2iX) 



r(i -fc + zA) 



dX, 
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where k < ^. Using the Whittaker function transform we see that the formulas in Theorem 111 101 
are a generahzation of Koornwinder's formula |13| (5.14)], stating that Laguerre polynomials are 
mapped onto continuous dual Hahn polynomials by the Whittaker function transform. 

(ii) The product formula in Theorem 13.101 has a similar structure as the discontinuous integral 
for Bessel functions of Weber and Schafheitlin, see [IHl §13.4]. 

(iii) Theorem 13. 101 can be obtained as a limit case of a bilinear summation formula for Meixner- 
Pollaczek polynomials, see [HI Rem.3.2(ii)]. 
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